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Abstract. A measurement of the forward-backward charge asymmetry in top-antitop (tt) pair production in proton- antiproton 
(pp) collisions is presented. The asymmetry is measured for different jet multiplicities in the lepton+jets final state on 0.9 fb _1 
of data collected by the D0 experiment at the Fermilab Tevatron Collider. The result is sensitive to new physics, which is 
demonstrated by setting an upper limit on tt production via heavy neutral gauge bosons (Z'). 
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INTRODUCTION 

At the Tevatron, top quarks are mainly produced in 
top-antitop-quark (tt) pairs. Since the initial proton- 
antiproton (pp) state is not a charge conjugation (C) 
eigenstate, the final state is not expected to be symmetric 
under C either. 

The charge asymmetry is therefore not expected to 
vanish and can be translated to a forward-backward 
asymmetry. Because of limited statistics, instead of a dif- 
ferential measurement this analysis [ 1 ] measures the in- 
tegrated forward-backward asymmetry: 



Nf-N b 
Nf+Nb' 



(1) 



where Nf (N/,) is the number of forward (backward) 
events, meaning events, where the rapidity difference 
&y — yt — yi between top quark and antitop quark is 
positive (negative). 

At leading-order (LO) tt production is charge symmet- 
ric. However, starting at next-to-leading-order (NLO), an 
asymmetry arises from the inferences between produc- 
tion diagrams. A total asymmetry of 4-5% is expected 
10] . The contributions to A «, differ depending on which 
processes interfere. They are positive (6.4% NLO) for 
2^2 processes while they are negative (-(7-8)% NLO, 
-(3-0)% NNLO) for 2^3 processes, which result in at 
least one more jet in the final state iH3l. By measuring 
the asymmetry for different jet multiplicities, this analy- 
sis attempts to verify this behavior. 



MEASUREMENT OF THE 
FORWARD-BACKWARD CHARGE 
ASYMMETRY 

The presented analysis uses approximately 900 pb~' of 
data collected by the D0 experiment [5] at Fermilab and 
proceeds with the following steps. First a sample en- 
riched in tt events is selected, then each event is fully re- 
constructed using a kinematic fitter. After that, the sam- 
ple is fitted for sample composition and asymmetry si- 
multaneously. To suppress model-dependent effects, the 
result is not corrected for acceptance, instead a simple 
description of the accepted phase space is given. Like- 
wise, it is not corrected for reconstruction effects, but a 
dilution is given to be used with any model. 



Event Selection 

Each the top and the antitop quark decays to a W bo- 
son and a b quark. The subsequent decay of the W boson 
determines the event topology. For this analysis, we con- 
sider the decay channel with one W boson decaying to 
hadrons and the other one to leptons. The resulting final 
state consists of one high-p^ lepton, missing transverse 
energy tyj from the neutrino, two jets originating from a 
long-lived b hadron (b jet), and two light quark jets. This 
will be referred to as the lepton+jets channel. 

Based on the final state, events with >4 reconstructed 
jets are selected. The jets have to have transverse mo- 
mentum pr > 20 GeV (35 GeV for the leading jet) and 
pseudorapidity tj | < 2.5. At least one of the jets has to 
be identified as a b jet by a neural network tagging algo- 
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FIGURE 1. Fitted likelihood discriminant for events with >4 jets. Events reconstructed (a) as forward (Ay > 0) and (b) as 
backward (Ay < 0). 



rithm [6]. We also require exactly one isolated electron 
with pr > 15 GeV and \rj\ < 1.1 or one isolated muon 
with p T > 18 GeV and |tj| < 2.0. The event has to have 
a primary vertex, which fulfills certain quality require- 
ments, and $t > 15 GeV. 

After the selection, each event is reconstructed to a it 
hypothesis by a kinematic fitter J7|]. It varies the four- 
momenta of each object within their resolutions and 
minimizes a % 2 statistic under the constraint that the 
W boson masses are exactly 80.4 GeV and the top quark 
masses exactly 170 GeV. The ^-tagging information is 
used to reduce the jet combinatorics. 

Although it is also possible to determine the aysm- 
metry using only the lepton charge and direction, recon- 
structing the whole event doubles the sensitivity. 



Sample Composition and Asymmetry 

The largest background processes to 1 1 production are 
W+jets and multijet production. The size and asymme- 
try of the multijet background is estimated from data, 
using a sample with non-isolated leptons. To be able to 
discriminate from the W+jets background, a likelihood 
discriminant Jzf is defined, which is based on variables 
which are well-modeled by Monte Carlo, provide good 
separation between it and W+jets, and do not bias the 
rapidity difference |Ay|. The asymmetry of the W+jets 
background is suppressed by the kinematic fit, and the 
remaining asymmetry is estimated by simulation. The 
other small sources of background are taken into account 
in the systematic uncertainties. The likelihood distribu- 
tions are fitted simultaneously for sample composition 
and sign of Ay using a maximum likelihood method. Fig- 
ureQ]shows the outcome of the fit for events with >4 jets 
in data and Monte Carlo. 

The resulting asymmetry is Aft, = [12 ± 8(stat.) ± 
l(sys.)]% for events with > 4 jets. Table Q] shows the 



asymmetry as it was measured on data and as predicted 
by the generator MC@NLO [8] for different jet multi- 
plicities. 

TABLE 1. Asymmetry in percent measured on data and 
as predicted by MC@NLO for different jet multiplicities. 



N 



jet 



D0 data 



MC@NLO 



>4 12± 8(stat)±l(sys) 0.8±0.2(stat)±1.0(acc) 



= 4 19± 9(stat)±2(sys) 
>5 -16^(stat)±3(sys) 



2.3±0.2(stat)±1.0(acc) 
-4.9±0.4(stat)±1.0(acc) 



Acceptance and Reconstruction Effects 

The integrated asymmetry strongly depends on the 
region of phase space probed. The largest effect stems 
from requiring >4 jets above a certain pj threshold, 
which is strongly correlated to the 4th highest particle 
jet pr- Figure [2] shows how the generated asymmetry 
depends on that variable. 
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FIGURE 2. The asymmetry as predicted by MC@NLO as a 
function of the 4th highest particle pj. 

Since the exact structure of the asymmetry in the 
Standard Model is unknown, we choose not to correct 
for the acceptance. Instead, the analysis is designed so 



that the selection can be described by simple particle- 
level cuts. This description has been shown to yield an 
accuracy of 2% (absolute). 

Misreconstruction of the sign of Ay dilutes the ob- 
served asymmetry. It can occur by misidentifying the 
lepton charge, a negligible effect, and by misreconstruct- 
ing the event geometry. If p is the probability to measure 
the sign of Ay correctly, the fraction of visible asymme- 
try, which is called the dilution factor, is D = 2p — 1. 
Since corrections for the reconstruction effects are highly 
model-dependent, the parametrization of the dilution for 
different jet multiplicities is provided instead. It is shown 
for events with >4 jets in Figure[3] 
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FIGURE 3. The dilution and its uncertainty as a function of 
generated |Ay| for //events with >4jets. 

To compare a predicted asymmetry from any model 
with the observed one, the generated A« has to be folded 
with the acceptance and reconstruction effects. This was 
done with the MC@NLO prediction to obtain the num- 
bers in Table Q] 

SENSITIVITY TO NEW PHYSICS 

Resonance production of tt can change the forward- 
backward charge asymmetry. Axigluons are predicted to 
give negative asymmetries [9]. Production via heavy neu- 
tral gauge bosons (Z'), however, can give large positive 
contributions to the asymmetry, and this example was 
studied in this work. By measuring the asymmetry one 
is in principle sensitive to narrow and wide resonances. 
The distribution for A«,, as a function of the fraction / 
of ff events produced via a leptophobic Z' with Z-like 
coupling to quarks, was predicted for several Z' masses 
using ensembles of simulated datasets. Employing the 
Feldman-Cousins method [10] limits at 95% confidence 
are derived and displayed in Figure|4] 
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FIGURE 4. 95% C.L. limits on the fraction of tt produced 
via a Z' resonance as a function of the Z' mass. The dashed 
curve indicates the expected limit with the shaded bands at 
one and two standard deviations. The solid curve represents the 
observed limit. 



the MC@NLO predictions, and the expected behavior 
for events with = 4 jets and > 5 jets is indicated by the 
data. Parametrizations of the acceptance and dilution al- 
low comparisons with any model. Limits on tt produc- 
tion via heavy neutral gauge bosons can be determined. 
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SUMMARY 



A first measurement of the integrated forward-backward 
charge asymmetry in top quark pair production was pre- 
sented. The measured asymmetries are consistent with 



